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Biomass energy systems can be employed to meet the requirements of distributed energy systems in 
rural as well as urban contexts, whether this is an electrification or a microgeneration project. This work 
is focused on a mathematical programming approach applied to bio-based supply chains that use locally 
available biomass at or near the point of use in order to produce electricity or other bioproduct. The 
problem of designing and planning a regional biomass supply chain is formulated as a MO-MILP (multi¬ 
objective mixed integer linear program), which takes into account three main objectives: economic, 
environmental and social criteria. The model supports decision-making about location and capacity of 
technologies, connectivity between the supply entities, biomass storage periods, matter transportation 
and biomass utilisation. The advantages of this approach are highlighted by solving a case study of 
a specific district in Ghana. The aim is to determine the most suitable biomass and electricity network 
among the different communities. The technology considered to transform the biomass into electricity is 
gasification combined with a gas engine. 

© 2012 Elsevier Ltd. All rights reserved. 


1. Introduction 

The development of a successful bioenergy sector in both 
industrialised and developing countries will contribute to a long¬ 
term diversity, security and self-sufficiency of energy supply [1], 
Current challenges emerging in the worldwide energy sector are 
consequence of diminution of natural sources, climate change and 
technology development. By adequately exploiting the resources 
and the available technologies, bioenergy is reckoned as (i) one of 
the most appreciated options to mitigate greenhouse gases emis¬ 
sions by replacing fossil fuels, (ii) an important resource in rural 
areas, and (iii) as a sustainable source of energy in vehicles and in 
electric power generation [2], 

Biomass markets are changing from exclusively local to inter¬ 
national markets. In this context, some regions are foreseen to 
become net exporters (e.g. Eastern Europe, Sub-Saharan Africa) and 
some other are potential importers (e.g. several members of the 
European Union, Japan) [3], In order to create a stable market with 
biomass being a commodity, efficient networks to sustainable 
supply the amounts of biomass required, cost effective technologies 
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to convert biomass, and improved distribution infrastructures to 
deliver the final product (i.e., energy or fuel) are to be developed 
[4], Even more, the biomass market also implies a social dimension 
which needs appropriate regulations and policies to motivate its 
implementation, as well as new behaviour patterns [5]. To face 
these challenges, biomass sector practitioners should have avail¬ 
able tools to guide the energy sector transformation in industri¬ 
alised and developing countries. Additionally, such tools should 
integrate sustainability criteria to support the decision-making 
process. In this context, a SC (supply chain) modelling approach 
provides the framework to adequately exploit biomass for energy 
purposes and to properly assess the trade-off among different 
contrasting metrics and different biomass sources usages. 

SCM (Supply chain management) deals with the coordination of 
materials, cash and information flows along the different processes 
in order to deliver efficiently goods or services. The processes and 
SC components can be geographically distributed. This is especially 
important for biomass to energy projects, that can be highly 
geographically dependent and their profitability can be strongly 
influenced by the location of the different processes and biomass 
sources. Indeed, it is widely recognised that logistic variables are of 
special complexity in this sort of projects [6], Commonly, biomass 
production and transportation account for a significant fraction of 
the whole bioenergy SC cost [7], SC modelling provides a tool 
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capable of evaluating the trade-offs between the different potential 
locations for processing and consumption sites, and the different 
feedstock sources and their specific properties, i.e. moisture 
content, energy and bulk densities, which as previously stated, is 
a requisite to develop efficient bioenergy networks. 

With regard to the performance evaluation, the implementation 
of a framework that considers the different realms of sustainability, 
is a challenge. This calls for a multiple-objective optimisation 
approach which should include economic, environmental and social 
criteria. Among them, the social or societal criterion is the less 
studied. Many rural projects traduce the social factor to exclusively 
job generation [8,9], The work by [ 10] defines a social function based 
on employment and liveability. Nidumolu et al. [11] include the 
maximisation of food production as a surrogate criterion. The papers 
by [12] and [13] describe a methodology to determine the societal 
LCA (life cycle assessment), aiming to define a standard parameter. 
Job creation, local/national recruitment, generation of employment 
and conditions of work are some of the issues to be evaluated in this 
societal LCA. Note that social concerns in industrialised and devel¬ 
oping countries can involve completely different aspects, ranging 
from security patterns to satisfaction of basic needs [14], 

1.1. Distributed energy systems 

Biomass as energy source has intrinsic characteristics that 
derive into technological limitations. Consequently, 100% biomass 
to energy projects employ small scale conversion systems which 
are generally placed close to the biomass generation source and/or 
to the demand points. This is to avoid logistic costs associated 
with the typical low energy densities and high volumes of biomass 
[6,15,16], Bouffard and Kirschen [17] point out that the electricity 
generation sector will have to combine the best attributes of both 
paradigms, decentralised biomass conversion systems and large 
power plants with biomass co-used with fossil fuels. 

Centralised energy systems consist in large power plants that 
export electricity to the grid and transport the raw material or 
energy source to the plant, while decentralised or DES (distributed 
energy systems) entail localised electricity generation near the 
demand points and near the biomass generation places. The 
conditions of the region makes either option advisable. No agree¬ 
ment exists in the literature about the definition of distributed 
generation, nevertheless it is usually perceived as small scale 
electricity generation [18], The reviews of [19] and [18] point out 
that the term can be referred to (i) stand alone applications, (ii) 
standby sources that supply power during grid outages, (iii) co¬ 
generation or waste heat recovery installations with power injec¬ 
tion to the grid, (iv) DES that support the grid by decreasing power 
losses and improving the system voltage profile, and (v) energy 
systems connected directly to the grid that sell the electricity 
produced. This work considers DES as stand alone applications. 
Decentralised energy projects in development contexts are char¬ 
acterised by a lack of an integrated approach in rural electrification 
planning [9], This has motivated the use of operations research to 
elucidate sustainable renewable energy approaches. For instance, 
[9] evaluates different technological options through goal 
programming to assess a qualitative response with respect to 
electricity generation cost, employment generation, avoided 
emissions and land use. Silva and Nakata [20] use LP (linear 
programming) to deal with the energy planning, finding financially 
sustainable solutions that appropiately connect sources and clients, 
and considers reliable income levels and energy demands. The 
subsidised costs are minimised. The paper suggests that the 
proliferation of renewable sources should be derived from envi¬ 
ronmental policies. Kanasse-Patil et al. [21] also employ an LP 
formulation to assure a reliable integrated renewable energy 


system by evaluating the cost of the energy, the costumer inter¬ 
ruption costs, and the expected energy not supplied. The paper by 
[22] takes into account a large number of state-of-the-art evalu¬ 
ating parameters which are used for decentralised energy planning. 
Goal programming is employed using cost, system efficiency, 
petroleum products usage, locally available resources, employment 
generation, emissions and reliability on renewable energy systems 
as objective functions. The results show that, in the considered 
context, biomass-based systems have the potential to meet rural 
needs. Furthermore, they offer reliability, promote local participa¬ 
tion, local control and create skills. The works by [23] and [24] 
develop, describe and use a multi-criteria decision tool called 
SURE, that is set to choose the appropriate energy mix of technol¬ 
ogies to match the energy demand of a rural area. The tool 
combines quantitative and qualitative parameters and allows for 
changes on the priorities according to the decision-maker criteria, ft 
is pointed out that technology assessment methods should be 
further developed to formulate more appropriated implementation 
strategies. The paper by [25] is an example of a modular renewable 
energy source implementation. It takes into account the results of 
a MILP approach that optimise the location of wind generators and 
selects the appropriate microgrid extension by minimising the 
investment costs. The approach for electricity generation through 
forest, fossil and wood-waste fuels described in [26] proposes 
a dynamic model to take operational decisions related to the 
biomass seasonality and availability. The interrupted availability of 
agricultural products has driven to an uncertainty treatment with 
a process graphs method [27], 

The present work is specially focused on regional biomass SC. 
There are different works dealing with this topic in the literature, for 
instance [28] and [29], In these cases, clusters and zones of 
competition and co-operation are determined using spatial plan¬ 
ning and LP. Fisher and Pigneri [30] identify potential feedstock 
supply sites and specific areas of electricity demand with geographic 
matching, evaluating the levelised cost of the energy. It is pointed 
out that techno-economic and social-institutional factors are critical 
in regional problems. Silva and Nakata [31 ] introduce the concept of 
regional disparity to distinguish among three possible situations: 
urban, rural and remote areas. Urban areas are defined as large 
cities: rural and remote areas concern small towns and dispersed 
population, whether or not connection to the grid is available. 

1.2. Small scale gasification 

Biomass can be transformed into a gas, by means of digestion or 
gasification: into a liquid through fast pyrolysis, or directly com¬ 
busted as a solid. In this paper, gasification is the chosen techno¬ 
logical option to provide a combustible gas. According to [2], large 
gasification systems are from 10 MW®, while small gasification 
systems cover the range from less than 100 kW t h up to a few MW*. 
In terms of electricity and in accordance with [32], small scale 
gasification plants enclose plants up to 200 kW e . 

Small scale gasification systems are used to meet the specific 
requirements of DES using the available biomass, or pre-treated 
biomass, such as pellets, bought from an established market. 
Gasification systems have to overcome several issues before being 
a well extended technology. Lessons learned from real experiencies 
reflect the need of improving the technology and enhancing 
economic efficiencies, as well as finding appropriate application 
domains. There is a lack of well trained technicians in the field and 
of monitoring and evaluation programs of biomass gasification at 
small scale plants. There is also a lack of focus and proper impetus in 
such a projects. Current challenges on small scale gasification plants 
are raw material characterisation and homogenisation, improved 
gasification operation results (concerning tars and producer gas 
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composition), more efficient cleaning units and an overall optimal 
performance [33], Nevertheless [34], describes how small scale 
biomass gasification has been successful in a few niche markets. 

The case study of this paper exemplifies a rural electrification 
problem. It tackles the strategic-tactical problem associated to the 
optimal design and operation of a BSC (bio-based supply chain) 
using a multi-objective mixed integer linear program (MO-MILP). 
The MILP approach captures the relevant fixed and variable oper¬ 
ating costs for each facility and each major product [35] considering 
equipment allocation. The SC problem defined includes a specific 
region, understanding region as a delimited area where self- 
sufficient resources and clients are embraced. The gasification 
system is linked to an internal combustion engine (G-ICE) to 
produce electricity. This approach differs from the previous works 
in the consideration of biomass properties changes along the SC, as 
well as in the evaluation of economic, environmental and social 
trade-offs, aiming to propose sustainable solutions in time, far from 
the usual cooperative solutions proposed for rural electrification 
projects in developing countries. 

2. The bio-based supply chain 

This paper proposes a BSC based on two blocks: the biomass and 
the energy blocks. The first one comprises sourcing, biomass pre¬ 
treatment and raw matter transportation, while the second one 
considers product generation and distribution. The sequence of 
pre-treatment, storage and transportation may change depending 
on the biomass type, the specific case study conditions and the 
chosen supply strategy. 

Special features that distinguish a BSC from a conventional SC 
problem are: the use of multiple biomass sources from different 
origins that are geographically distributed, and the subsequent pre¬ 
treatment required to homogenise the material in mass and energy 
terms. These features imply the combination of different MC 
(moisture contents), DM (dry matters), LHV (lower heating values) 
and BD (bulk densities). Biomass, with high MC, low BD, low LHV 
and fibrous nature, may lead to biomass pre-treatment so as to 
optimise its transport, handling and processing. Biomass properties 
can change along the SC, even if it is not the desired effect. The 
abovementioned biomass attributes are characterised in the SC 
model presented herein. 

The major steps in the BSC are described as follows: 

• Biomass growing, harvesting and collecting involve biomass 
production, by recovering biomass waste or using energy crops. 
Processes included here are drying, i.e. natural drying in the 
land field, baling or chipping. The two latter processes are used 
to diminish the volume and reduce the risk of deterioration 
[36]. Resource seasonality determines the harvesting or 
collection period. Different seasonal sources mix and storage 
could mitigate the impact on supply continuity. 

• Biomass pre-treatment includes all the necessary steps to 
produce an upgraded fuel. Such a fuel is homogeneous, has no 
impurities and is denser in terms of mass and energy. Pre¬ 
treatment allows costs reduction in biomass treatment, 
storage and transportation. Briquetting, pelletisation, torre- 
faction, pyrolysis, and combinations of them, are pre-treatment 
techniques [3], Research in biomass pre-treatment is crucial for 
the large scale biomass market. In this line, [7] contemplate 
forest wood residues to produce torrified wood that supplies 
gasification units to produce electricity. Magalhaes et al. [37] 
evaluate prices, yields, distances, investment and operation 
costs of torrefaction and fast pyrolysis in a real case study that 
uses imported biomass. The paper by [38] evaluates the SC of 
bioslurry, a mixture of bio-oil and biochar, through the trade¬ 


off between the investment cost of distributed pyrolysis units 
and the reduced operation costs due to a denser biomass. 
Storage can be considered throughout the BSC. Biomass can be 
stored whether before the transportation stage, at the biomass 
origin, in an intermediate step or at the power station site. This 
is crucial to create backups against seasonality. A drying 
phenomenon takes place here, even some dry matter loss may 
occur [39], Costs depend on the location and type of storage, i.e. 
open air, roof covered, air fan, indoor storage. Their selection 
depends on climate conditions, shape and volume of biomass 
and time of storage [16], The work by [39] emphasizes the 
multi-biomass seasonal availability combined with the 
biomass storage problem. 

Transportation can be expensive due to the biomass low energy 
density. Costs depend on distance, speed, tortuosity, hauliers 
capacity and amount of biomass to be transported. Operating 
costs such as driver remuneration or fuel consumption, as well 
as social and environmental impacts should be evaluated [16]. 
Forsberg [40] analyses the transportation problem in terms of 
environmental viability and ecological sustainability. It 
considers the utilisation of trailer, truck, train, ship or pipeline, 
thus implying different biomass states and different distances 
to cover. Optimisation of this activity implies to find the 
appropriate match between the amount of disposable biomass 
near the plant and the plant size. Yu et al. [41 ] evaluate a malle 
residues SC focused on transportation costs. It results that on- 
farm transportation for central biomass gathering, can be 
more expensive than biomass road transportation due to the 
elevated costs derived from the feedstock collection area. The 
papers by [42] and [43] compare the transportation of bio-oil by 
pipeline or by truck through an LCA. The results depend on the 
distance to be covered and the source of energy consumed to 
pump the bio-oil or to run the trucks. Depending on the scale, 
transportation could represent the bottleneck of the BSC. 
Biomass treatment refers to the biomass-processing plant 
to produce the desired product, i.e. biofuels, bioproducts, 
heat, cool or electricity. Multi-products can be also considered, 
as in [44], 


3. Problem statement 

The design-planning of a BSC may require a highly application- 
specific or tailor-made approach [45], The optimal SC specifically 
fulfills the needs of the target market; thus it is not possible to 
scale-up the proposed solution to address a different problem. 

3.1. Inputs 

1 Process data 

• A set of materials: raw matter, intermediates and final 
products. 

• A set of demands. 

• A proposal of biomass pre-treatments and treatment tech¬ 
nologies characterised by their efficiency. Consequently, 
a set of biomass materials as raw, pre-processed and pro¬ 
cessed matter, depicted by their main properties, are 
determined accordingly. 

• A set of providers, intermediates and consumers locations. 

• A time horizon. 

2 Economic data 

• An interest rate. 

• Investment, fixed and variable costs associated with all the 
technological options involved. 

• Products and consumables prices. 
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3 Environmental data 

• Raw material production environmental interventions. 

• Each process environmental interventions. 

• Transportation environmental interventions. 

4 Social data 

• Number of communities and number of potential processes 
to be installed in each community. 


3.2. Outputs 

1. The most suitable pre-processing units, with their corre¬ 
sponding capacity and location. 

2. Connections between providers, intermediates and consumers. 

3. Biomass storage periods. 

4. Matter transportation flows. 

5. Biomass utilisation. 

6. The detail of the economic parameters and environmental 
impacts for each SC echelon. 

The chosen objective functions should represent the decision¬ 
maker interests. The metrics considered in this work include the 
three branches of sustainability. 

• Net present value (NPV) as the economic metric. 

• An LCA is performed for the environmental issues. 

• The creation of jobs in the widest range of communities as the 
social indicator. 

4. Mathematical model 

A general outline of the BSC is shown in Fig. 1. It contemplates 
the activities described in Section 2: sourcing, pre-treatment, 
product generation and product distribution. 

The BSC is defined as a number of potential locations where 
processing sites or distribution centres, or both of them can be 
installed. Suppliers are at fixed locations where biomass is avail¬ 
able. The final product can be produced at several sites. The prop¬ 
erties of the raw biomass may be enhanced by means of the pre¬ 
treatment units so as to allow intermediates to meet the charac¬ 
teristics required by subsequent steps in the SC. Pre-treatments 
may be also convenient for reducing transportation costs. The 
production capacity of each processing site is modelled by relating 
the nominal production rate per activity to the availability of the 
equipment per year (i.e. the number of working hours per year and 


equipment). Storage and transportation activities are modelled by 
considering upper and lower bounds based on their corresponding 
biomass handling capacity. Distribution centres can be supplied 
from more than one pre-treatment plant or supplier. Materials flow 
between facilities may appear if selecting such a flow improves the 
performance of the SC. The market demand of electricity can be 
satisfied by more than one site. 

All the abovementioned decisions will be taken by considering 
the optimisation of an economic (NPV), an environmental impact 
(Impact 2002+) and a social metric (SoC j along a pre-defined 
planning horizon. 

Next, the mathematical formulation is described. The resulting 
model is solved by using an MO-MILP technique, specifically the 
e-constraint method. Multi-objective optimisation allows assessing 
the trade-off among the economic, the environmental impact and 
the social indicator. The model variables and constraints can be 
categorised into four groups: (i) process operations constraints 
given by the design-planning sub-model, (ii) the economic metric 
formulation, (iii) the environmental sub-model and (iv) the surro¬ 
gate social metric. 

4.1. Design-planning model 

The design-planning model selected to deal with the BSC 
network is adapted from the work of [46], This model translates the 
State-Task-Network (STN) formulation [47], which is a widely 
known approach for scheduling, to the SC context. This SC design¬ 
planning formulation can collect all the SC nodes activities infor¬ 
mation through a single variable set. On the one hand, this feature 
eases the economic and environmental metrics formulation. On the 
other hand, it facilitates the consideration of pre-treatment activi¬ 
ties and their outputs. The SC material balances can be modelled by 
means of a single equation set for all materials and echelons. This is 
possible since manufacturing nodes and distribution centres, as 
well as final products, raw material and intermediates are treated 
indistinctively. Therefore, the most relevant variable of the model is 
Pyff't [kg]- which represents the magnitude of a particular task i, 
performed using technology j during period t, whose origin is 
location / and destination is location /'. In the case of production 
activities, they must receive and deliver material within the same 
site [Pijfft). In contrast, in a distribution activity, facilities / and /' 
must be different. This mathematical formulation assumes that an 
activity consumes and produces certain materials with determined 
properties and can be performed in different equipment. Using the 
activities as the core of the formulation rather than products- 





Fig. 1. General scheme 
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materials renders a flexible formulation which can be easily 
extended to deal with different case studies. The equations 
comprising this formulation are described in the following 
paragraphs. 

Mass balance must be satisfied at each node of the network. The 
expression for the mass balance for each type of material s (raw 
material, pre-processed biomass, final product) processed at each 
potential site / in every time period t is presented in Eq. (1). 
Parameter a s y [adim.] is defined as the mass fraction of material s 
that is produced by task i using technology j. T s is a set that refers to 
tasks that produce s, while a s y and T s sets, are associated with tasks 
which consume s. 

= EE E a sij P ijf'ft 

f ' d) 

-EE E a sijPimVs,f,t 
r 

The change in inventory of material s for consecutive planning 
periods and for each location is given by the difference between the 
amount of material produced by those tasks belonging to set T s and 
the amount consumed by those tasks included in set T s . The model 
assumes that process parameters such as conversions, separation 
factors or temperatures, are fixed for each activity in order to 
enforce the linearity of the problem. In this sense, the parameters 
a S jj [adim.] and a s y [adim.] give the “recipe" for a specific activity. 
Nevertheless, there are activities (7) for which it is desirable to let 
the model specify the mixture of inputs in order to achieve a given 
value of a specific biomass property, for instance, a specific mois¬ 
ture content. For such activities, the combination of feedstocks and, 
therefore, the proportion of each feedstock is variable. In order to 
take into account such activities, the mass balance is modified as 
shown in Eq. (2). Note that Eq. (1) is a particular case of Eq. (2) 

Ssfi-Stft 1 =EE E a sij P iff'ft “EE E 

f ' ieTs J-(/«n]/0 f 

+ E E Pv m 

ie(%ni) jE fanjj,) 

-EE Pv m (2) 

fe(T,n/) ie ( 7i nJ r ) 

With regard to the variables Pyff, t [kg] and Pv si jf, t [kg], notice that 
the former is used in the mass balance (Eq. (2)) coupled with the 
parameter a si j or a s y which specify a fixed proportion of material 
produced or consumed for a task i. On the other hand, Pv^ is 
modelling flexible tasks which allow the proportion of the material 
produced or consumed to vary so as to provide more alternatives 
for biomass mixing to the model. For this reason, Pv si jft is not 
multipled by such parameters in Eq. (2). 

For the activities for which biomass properties are fixed, the 
energy balance is satisfied directly by the definition of the streams. 
However, it is necessary to verify that the energy balance is satisfied 
for the flexible activities (7). The energy balance for such activities is 
represented by Eq. (3). Here, HV S [MJ/kg] is the heating value of 
material s. Notice that each different type of biomass has a different 
heating value. A specific activity changes the heating value of the 
output stream if (i) it is a pre-treatment task that modifies explicitly 
the calorific value of the biomass, or (ii) it is a task whose main 
objective is the change of shape, but it is receiving a mixture of 
biomasses as input. 

E WsPvsijft = E ™s p v m Vt e 7, j, f,t (3) 

KT, seT, 


In case the flexible activities must accomplish an MC for the 
input stream, constraint (4) must be satisfied. The parameters 
Water s [adim.] and Water]™* [adim.] represent the MC for material 
s, and the maximum MC allowed for task i performed in equipment 

j, respectively. 

E Water s Pv sljft < Waterfi 10 * ]T Pv sijf , t V ie I, j, f,t (4) 

seS ‘ se 5,- 

Another important extension of this model is the consideration 
that storage is capable of changing biomass properties (e.g., mois¬ 
ture content, dry matter, heating value). In order to do so, storage 
should be considered as an actual activity. Let us define the subsets 
Jstor and S s tor which will represent the storage “equipment” and 
those materials that when kept in storage change their properties, 
respectively. Notice that the mass balance has been decomposed in 
this case so as to deal with the one period delay necessary for the 
properties change to occur. As expressed by Eqs. (5) and (6), 
a storage activity places inventory in the current period t and takes 
inventory from the previous period t - 1. 

% = EE E «$ p m v^s stor ,f,t (5) 

f iGT ’ jc (jiCijf, nw) 

% 1 = E E E ^si] p ijff't VScW,D t (6) 

f ' JE (jiCi) f n] swr ) 

We should point out that an important limitation of this 
modelling approach is that Eqs. (2)—(3) allow combining different 
materials in proportions which are defined by the optimisation just 
for the activity which performs the energy generation. The states 
feeding any other activity and their corresponding proportions and 
properties (e.g. MC, DM, BD, LHV) must be defined a priori. This 
implies that no combinations of different biomass are allowed 
before storage or such combinations must be also defined a priori. 

Eqs. (7) and (8) represent the temporal change in the equipment 
technology installed in the potential facility locations. We will 
considered economies of scale by using a piecewise linear 
approximation in K different intervals and a so-called SOS2 variable 
type (Zjjk). Such variables are positive and at most two consecutive 
variables are non-zero. EEjjj™ t [h] is the limit of capacity for interval 

k. Vjf t is a binary variable indicating whether or not capacity of 
technology j is expanded at site/in period t. This formulation will 
be recalled in the economic performance metric section for 
computing the investments associated with capacity expansions. 
Eq. (9) is used for the total capacity F]f t [h] bookkeeping taking into 
account the capacity augment during the planning period t (FEjg 
[h]) for equipment technology./ in facility/. This equation considers 
the case of the initial design of an SC (FEjf 0 — 0) as well as an SC 


retrofit scenario {FEjf 0 j= 0). 


EfotHS* = V /.M7M 

k 

(7) 

E^ = % v/.jWr, t 

k 

(8) 

% = %- i+ f % v /> H f , t 

O) 


Eq. (10) is used to ensure an utilisation greater than or equal to 
a minimum value established by the decision maker and that the 
utilized capacity is lower than or equal to the available one. 
Parameter [adim.] defines a minimum utilisation rate of tech¬ 
nology j in site / as a proportion of the total available capacity. 
Parameter dyjj, [h/kg] represents the resource utilisation factor. This 
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is the capacity utilisation rate, in terms of capacity units (e.g., 
machine-hours), of technology j by task i whose origin is location/ 
and destination location/'. 

1<EE e w p w't < Fin -1 v/, m /5 t (io) 

f fe/j 

The capacity is expressed as equipment j available time during 
one planning period, then 6^ represents the time required to 
perform task i in equipments per unit of produced material. Thus, 
once operation times are determined, this parameter can be readily 
approximated. 

Eq. (11) guarantees that the amount of raw biomass s purchased 
from site/at each time period t is lower than an upper bound given 
by physical availability A s f t [kg] (e.g., seasonality, crop/plantation 
yield in a specific region). 

E E EVr<4/ t VscRM, feSup, t (11) 

f ieT s M 

Eq. (12) establishes that flows of energy exist only if locations/' 
and/are interconnected. Zf,f is a binary variable which takes a value 
equal to one if/' and/are interconnected, 0 otherwise; while M 
represent a big positive number. By Eq. (13) sales of final product 
se FP carried out from facility location /' to market /e M are esti¬ 
mated. Eq. (14) is used to express that the demand can be partially 
satisfied due to biomass production or supplier capacity limita¬ 
tions: the sales of product s carried out in market/during the time 
period t could be less than or equal to the demand. 

Pijfft < MZf, f \/SeFP, ie(T s CTr), feMkt, t (12) 

Sales sf , ft = ]T P m VseFPJsMktf <£Mkt,t (13) 

i e (T s nrr) je (; i nj / ) 


]T Sales sf ,jt < Dem S ft VseFP, feMkt, t (14) 

f'£M 

Here, it is important to emphasize two aspects. One is that for this 
sort of networks the final product is the energy delivered to the 
different regions (MkT). Thus, for those tasks that carry out the energy 
generation their corresponding parameter a s ,j is determined as 
a function of the heating value of the input materials and the efficiency 
of the equipment. Finally, the energy consumed by the equipment 
along the SC should be discounted from the energy available to satisfy 
the demand once the network reaches “steady state”. 

4.2. Economic model 

The expressions required to compute the operating revenue, the 
operation costs, the total capital investment, and NPV are next 
described. 

The operating revenue is expressed in Eq. (15) as the product 
sales during period t. 

ESalest = E E E Sales sf , ft Price sft Vt (15) 

seFP feMkt f'<£(MktUSup) 

The operating costs include fixed and variable costs: Eq. (16) 
describes the total fixed costs of operating the SC network. FCFJj^ 
[US$/h] is the fixed unitary capacity cost of using technology j at 
site/. 

FCost t m J2 E FC %% Vt (16) 

f &(MktU Sup) j e ] f 


In turn, as variable costs, the cost of purchases from supplier e, 
includes raw material procurement, transport and production 
resources, as shown in Eq. (17). The purchases of raw materials 
(Purch^p [US$]) made to supplier e are evaluated in Eq. (18). We will 
assume a different supplier for each component of the variable 
costs. This assumption can be easily relaxed to account for the 
specifics of the problem being dealt with. The variable Xest [US$/kg] 
represents the cost associated with raw material s purchased to 
supplier e. Transportation and production variable costs are 
determined by Eqs. (19) and (20), respectively. p^, t [US$/kg] 
denotes the e provider unitary transportation cost associated with 
material distribution from location/to location/' during period t. 
T gfet [US$/kg] represents the unitary production cost associated to 
perform task i using technology/ whereas [US$/kg] represents 
the unitary inventory costs of material s storage at site /. The 
parameter Tj^J t entails similar assumptions to the ones considered 
with regard to a si j and a si j, namely, the amount of utilities and 
labour required by an activity are proportional to the amount of 
material processed. 


EPurchet = Purch™ + Purch£ t + Purch.P r t od Ve, t 

(17) 

Purch™ = E E E T, P m*est VecE r m. t 

sg RM feF e ie % M 

(18) 

Purch* = E E E J2 P m'tPeff't VeeEtf, t 

, e Fr je Ji nj e / f 

(19) 

Purchr=EE E P m^t 

f i£Tr j e (j irt j f ) 

(20) 

+ E E S sft T sfet VeeEprod, t 

s fe(SupUMkt) 


The total capital investment on fixed assets is calculated by 
means of Eqs. (21) and (22). These equations include the invest¬ 
ment made to expand the technology’s capacity j in facility site/in 
period t. The investment required to connect two different loca¬ 
tions/and/' by using a medium voltage network (Invest [US$/ 
km]) is just accounted in the first planning period. Recall that 
economies of scale for technologies capacity is considered using 
a piecewise linear approximation in K intervals. Here, Price 
[US$] is the investment for a capacity expansion equal to the limit of 
interval k (FEj™" [h]). 

FAssett = E E 'E Price ifk it %m 

f j k (21) 

+ E ln vest MV distancefj.Zjj, Vt = 0 

f f 

FAssett = E E E Price jfk it '-jfkt Vt>0 (22) 

/ j k 

Finally, Eq. (23) calculates the profit in period t, as operating 
revenues minus fixed and variable operating costs. The NPV is 
calculated as in Eq. (24). 

Profitt = ESalest - f FCost t + y EPurch et J Vt (23) 


NP V = vY mfitt ~ FAsset t ) 

P+ir/ / 


(24) 
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4.3. Environmental model 

The application of the life cycle assessment (LCA) methodology 
to the SC model requires four steps according to ISO documents 
(ISO14040): goal definition and scope, life cycle inventory (LCI), life 
cycle impact assessment (LCIA) and results interpretation. Envi¬ 
ronmental interventions are translated into metrics related to 
environmental impact as end-points or mid-points metrics by the 
usage of characterisation factors. Eq. (25) calculates IC a ft [points] 
which represents the mid-point environmental impact a associated 
with site/, as a consequence of carrying out activities in period t. In 
turn, \l/ijjj, a [points/kg] is the a characterisation factor of the envi¬ 
ronmental category impact for task i performed using technology j, 
receiving materials from node/and delivering them at node/'. 

% = EEEWp (25) 

Mr ^ f 

Analogously to a S y and a si j, the value of ^ijff'a is fixed and 
constant, since all environmental impacts are considered linearly 
proportional to the activity performed in the node (variable Pij/ft 
[kg]) [48], Environmental impacts associated with transportation, 
has as a common functional unit (FU), the amount of kg of trans¬ 
ported material over a given distance (kg-km). Consequently, the 
value of the mid-point environmental impact ^ijjp a associated with 
transport, is calculated as in Eq. (26), where 4yff'a [points/(kg km)] 
represents the a characterisation factor of the environmental 
category impact for the transportation of a mass unit of material 
over a unit of length. Note that the impact is assigned to the origin 
node. The environmental impacts associated with production (Eq. 
(25)) or transportation (Eq. (26)), can be computed by setting the 
indices summation over the corresponding tasks (ie Tr or ie NTr). A 
tortuosity factor (Tortuosity [adim.]) may be employed to correct 
the estimated distance between nodes. 

'Pijjfa = ^Jjadistanceff,Tortuosity VieTr, jeJi, a, f, f (26) 

Eq. (27) introduces DamCgft [points], which is a weighted sum of 
all mid-point environmental interventions. They are combined 
using g endpoint damage factors ‘Cag [adim.], normalised with 
NormFg factors. Moreover, Eq. (28) calculates g normalised endpoint 
damage along the SC (DamC^ c [points]). 

DamCgft = NormF g Z ag IC aft Vg, /, t (27) 


DamCg C = ^ DamC glt V S ( 28 ) 

/ t 

Eqs. (29) and (30) sum the endpoint environmental damages for 


each site/and for the whole SC, respectively. 

Impactj 002 = J2 E DamC gft v / 

(29) 

Impact™ 2 all = E E J2 DamC gft 

f s t 

(30) 


For further details regarding the operational and environmental 
formulation the interested reader is referred to [46] and [49], 


4.4. Social model 

The proposed approach consists in a criterion that counts the 
number of sites that have a treatment or pre-treatment system 
installed. The aim is to install as many as possible to promote 


working places in the widest range of communities or demand 
sites. Therefore, the social criterion SoC should be maximised (see 
Eq. (31 )). This criterion assigns a value of 1 to each unit installed per 
site / V is the binary variable that characterises the number of units 
installed per site. This is not applied for storage. 

SoC = EEE% *;./>* ( 31 ) 

j r t 

SoC can be used as objective function along with the NPV and 
Impact in the MO-MILP formulation. An e-constraint approach 
will be used to tackle the MO-MILP. 

The overall optimisation problem can be posed mathematically 
as follows: 

Mm | - NPV, Impact™ 2 a , -Soc} 
subject to Eqs. (1)—(31); 

Jf’e (0,1}; yeR + 

Where, denotes the binary variables set, while y corresponds to 
the continuous variable set. 

4.5. Multi-objective strategy 

The multiple criteria decision-making process aims at finding 
the most suitable solution for the decision-maker point of view. 
Decision problems usually present multiple and conflicting criteria 
to evaluate the available alternatives. It is then necessary to make 
compromises or trade-offs regarding the results of the different 
possible choices. Instead of giving one specific solution, the multi¬ 
objective optimisation provides the so-called PF (Pareto frontier). 
This frontier represents a set of solutions that are not non- 
dominated by any other feasible solution [50], This paper 
assumes that the decision-maker priorities are not known, there¬ 
fore, the analysis is focused on a general perspective based on the 
characteristics of the PF and the most important features of the 
scenarios that optimise each one of the selected objective 
functions. 

• Concerning individual objective functions optimisation, 
extreme scenarios are obtained and characterised by solutions 
derived from NPV maximisation, Impact 2002 + minimisation 
and SoC maximisation. As more than two objectives are 
considered, iso-lines are used to graphically represent the 
three parameters bidimensionally. The iso-lines represent 
a fixed value for the SoC range. This is possible due to its 
discrete nature. The range of values is determined using the SoC 
values obtained from the extreme scenarios. 

• Once the SoC range is established, for each of its discrete values, 
a bi-objective optimisation for the NPV and Impact 2002+ is 
carried out. Firstly, two extreme points are obtained by the 
single optimisation of NPV and Impact 2002 + considering 
a fixed value for the SoC. The intermediate PF points are then 
determined. Those Pareto solutions are found using 
e-constraint method, that utilises a set of single objective (NPV) 
optimisation cases, while the remaining objective (Impact 
2002+) is set as an upper limit constraint. In this specific case 
study, 10 intermediate points are determined. Each one of the 
12 solutions is a different network configuration. Please find 
more details about the e-constraint methodology in [51]. 

• One of the disadvantages of some multi-objective optimisation 
methodologies such as the e-constraint and the weighted sum 
is that they may provide dominated solutions, i.e., solutions for 
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Fig. 2. Location of Atebubu district, in Ghana. 


which other ones that improve an objective without worsening 
another one are available. A filtering step is applied to remove 
such solutions to construct the PF. The filter works by 
comparing every point with every other generated point. If 
a point turns out to be dominated by other point, it is elimi¬ 
nated [51]. The decision-maker will then select a point 
(network) belonging to the PF. 

5. Case study: A bio-based supply chain located in Ghana 
using gasification 

The generic BSC considered in the mathematical formulation 
(Fig. 1 ), is adapted to suit the intrinsic nature of the case study: 

• A regional configuration deals with short distances. Owing to 
this, chipping and drying are the selected pre-treatments. 
Advanced pre-treatment methods such as fast pyrolysis, tor- 
refaction or pelletisation are disregarded since the small 
biomass quantities required to supply the electricity demand 
does not allow to take profit from economies of scale. In 
addition to that, they are too complicated technologies for rural 
areas in developing countries. Biomass waste is pre-processed 
before gasification to obtain the requiered shape, LHV and MC. 
Storage serves as a biomass reservoir to leverage capacities 
from several sites and seasonality. 

• The region does not have access to the grid. Customers can be 
therefore connected to a specific-built LV (low voltage) or a MV 
(medium voltage) microgrid. The LV distribution line supplies 
a community and the MV line, of 36 kV, connects different 
communities. The investment cost is contemplated. 

• Each community represents one unique consumer. 

This case study contemplates a specific rural area of a devel¬ 
oping country, Ghana (Africa). Nine communities in Atebubu- 
Amantin district, in the Brong Ahafo Region (see Fig. 2) are of 
concern. The biggest community is Kumfia with 2834 people, and 
the smallest one is Nwunwom with 122 people. See Appendix A for 
further details. The raw material is cassava waste. Cassava is a well 
extented tropical crop in the country, mainly used to provide food 
(in the form of fufu or gari). Agriculture is the most important 
economic sector, with cassava, yam and maize as basic products. 

The selected communities form a region with a main charac¬ 
teristic: they are equipped with a MFP (multi-functional platform) 
that currently supplies the electricity needs; cell phone’s battery 
charging, water refrigeration, lighting, radio, TV, computer and 
maize mills and cassava graters. Those communities are therefore 


used to pay for the electricity service. Data for the BSC character¬ 
isation are provided by local entities. The Energy Center, in Kwame 
NKrumah University of Science and Technology (KNUST), and the 
NGO Kumasi Institute of Technology, Energy and Environment 
(KITE). The MFP’s project in Atebubu district has been executed in 
collaboration with KITE and the NGO Women and Children Support 
Organisation (WACSO). Prices for transportation and utilities val¬ 
orisation, as well as specific information about the cassava crop, 
have been obtained during an on-field travel. The boundaries 
considered are from cradle-to-gate. The currency is US$2010- The 
applied currency conversion is 1 Ghana cedi per US$0.65946 [52], 


5.1. Raw material 

Cassava is also called manioc, manihoc, yuca, mandioca, aipin, 
catelinha, macaxeira and tapioca. The cassava waste is pilled in the 
farms and used for fertilisation or to feed pigs or poultry. Leaves from 
cassava crops are used for food (soups). The stalk is used to replant 
the crop. The underground cassava plant can be divided into 
rhizome and root. In turn, the root has three main components: pulp, 
starch and peels. $tarch is the main product and it is used for food. 
The most extended uses for cassava residues is ethanol production 
(pulp), biogas and activated carbon production (peels) [53,54], 

Cassava rhizome is the biomass waste considered in this study. 
The cassava is planted once a year, in April, during the rainy season 1 
and does not need any special care. It takes from 3 to 6 months to 
grow, depending on the plant type. The harvesting period is 
assumed to be from June to October. A 66.5% of the tubercle is 
cassava rhizome [55], As the produced wastes have no current 
alternative use, the cost of acquisition is considered negligible. 
Table 1 summarises the cassava main properties. The last column 
sums up the total amount of cassava waste produced by the 9 
communities. This value has been calculated by considering that 
from the total amount produced in Atebubu district, only a 20% can 
be taken for electricity purposes. According to [56], 264649t of 
cassava has been produced during 2009 in the whole Atebubu 
district, representing 17667.6 GJ/y. 

Fig. 3 represents the production of cassava waste per commu¬ 
nity in GJ/month. The agricultural activity is supposed evenly 
divided among the communities according to their population. 
Cassava waste disposal takes place at the same generation place. 
5ee also in Fig. 3 the aggregated thermal demand estimated for the 


1 Conversation with cassava farmers from Seneso. 
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Table 1 

Feedstock properties. 

Biomass Cost LHV„ MC Seasonality Yearly 

_($/t) (MJ/kg) {% wt)_available (t) 

Cassava waste 0 10.61 42.50 June-October 1666.13 


nine communities. These numbers are derived from the electricity 
demand and the gasification efficiency. 

5.2. Technologies used and electricity demand 

Inlet biomass to the gasifier must be chipped and have 20% or 
less of MC. It is supposed that the only possible biomass storage is 
carried out before chipping and gasification using on-field storage 
which is the cheapest and simplest option. Fig. 4 shows the layout 
of the different pre-treatment options applied to the biomass (BM). 
MC, DM, shape and LHV change along the network. 

It is assumed that chipper and dryer work an average of 8 h/d. 
The gasifier is assumed to work 16 h/d. The project has a lifespan of 
10 years [57], Estimation of chipper and dryer available range of 
capacity takes into account the processed maize in an MFP during 
one day and the total production of cassava estimated per year. The 
investment and O&M costs are taken from the literature. Economies 
of scale are considered. The percentages and proportions used to 
linearly model the activities in the mathematical formulation and 
to calculate the biomass states are summarised below. 

• Biomass waste generation. The cassava is harvested and the 
rhizome is separated from the rest of the plant. 

• On-site pre-treatment. Biomass is stored in biomass generation 
sites. Therefore, open air storage, in open air piles, takes place 
here. 

• Drying. Passive drying is a consequence of the open air storage. 
According to [58], storage at the roadside is free and has 
negligible O&M cost. Raw biomass with an MC higher than 20%, 
has a DM loss per month of 3%, and an MC loss per month of 2%. 
12 months is the maximum storage period. Active drying takes 
place into a rotatory drum. This unit is used to decrease the 
inlet MC to the required 20%. Biomass changes its MC and its 
LHV according to the tonnes of water evaporated. This unit has 
an energy efficiency of 99%. It consumes diesel as utility. The 
available capacities are between 0.1 and 5 t/h. 

• Chipping. Chipping is placed after drying. It consumes elec¬ 
tricity from the G-ICE system. This unit has an energy efficiency 
of 96%. The available capacities are in the range of 0.1-5 t/h. 

• G-ICE system. The most important parameters in the biomass 
gasifier, which influence the producer gas composition, are the 
amount of inlet air (i.e., the equivalence ratio, ER) and the MC. 
See the main system’s parameters in Table 2; temperature of 
the producer gas at the gasifier outlet stream {T gaS if), producer 
gas composition and LHV, power and heat outputs, as well as 
the partial (cold gas efficiency, CGE) and overall (77) efficiencies 
of the system [33]. p is used in the mathematical model to 
describe this activity. The gasification units range between 5 
and 100 kW e . 

• Transportation. Solid biomass can be displaced from its point of 
generation to a storage place or to a pre-treatment node by 
tractors, within a capacity of 10t. 2 Lineal distances, with 
a tortuosity factor of 1.8, are considered and expressed in km. 

• LV and MV distribution grids. The LV distribution line has a 6% 
losses in energy terms. The MV distribution line takes into 



account losses estimated by [59], which are proportional to the 
power demand. See Appendix A for further details. 

The economic evaluation is based on estimations of investment 
and annual O&M costs. The total capital requirement is spent at the 
beginning of the project. Table 3 lists the parameters needed for 
costs estimation of pre-treatment units, G-ICE plant and trans¬ 
portation. The diesel price is $1133.31 /t. 3 The electricity price is 
$0.233/kWh, as the value that the communities are willing to pay. 
Its calculation takes into account the consumption of diesel and 
batteries, and their respective prices based on data from the Seneso 
community. 

The electricity demand has been estimated on the basis of 
references from previous experiences on rural electrification 
projects, in West Africa and South American communities. It takes 
into account residential, community and commercial electricity 
requirements. See additional details in Appendix A. The highest 
gross demand is 448.65 lcWh/d in Kumfia community, while the 
lowest is 21.17 kWh/d in Nwunwom community. Such figures take 
into consideration the LV microgrid losses. Fig. 5 depicts the nine 
communities represented by their relative energy demands (in 
blue), in a square grid that represents their relative distance in km. 
Black points mark potential locations for pre-treatment and treat¬ 
ment units; all the communities and four more intermediate sites 
are considered as potential locations. 

5.3. Environmental impact 

Impact 2002 + metric evaluates the environmental impact in 
points (pts). LCI values are retrieved from LCI database [62] using 
[63], and they are directly converted into Impact 2002+ mid-point 
indicators in the so-called LCIA step. The evaluated BSC echelons 
are biomass production (cassava rhizome obtaining) without 
transportation, transportation by tractors, all pre-treatment possi¬ 
bilities (chipper and dryer) and generation of electricity by means 
of biomass gasification.The impact of this last task has been adap¬ 
ted from a large scale gasification plant impact, an IGCC, consid¬ 
ering efficiency difference, in energy terms, among the two plants. 
For further details on LCI values, see [33], 

5.4. Model characteristics 

This case study includes 40 different biomass states (s). Set j 
includes 6 different equipment technologies in accordance with the 
different pre-treatment and treatment technologies and means of 
transportation. The set of activities i comprises 79 elements 


2 On-field 


i WACSO. 
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Fig. 4. Pre-treatment activities layout. 


depending on the couple biomass-processing and transportation 
activities. Set / considers sites location, being overall 31 of them; 
involving suppliers, pre-treatment, treatment, potential locations 
and markets sites. This case study considers biomass evaluated in 
mass basis for transportation economic accounting. The planning 
period t represents a month. The project is evaluated along 
a planning horizon of 10 years. The interest rate is 15% [64], 

The mathematical model has been written in GAMS and solved 
using CPLEX 11.0 on a PC Intel(R) Core(TM) i7-2620M CPU 2.70 GHz 
and 4.00 Gb of RAM. The optimisation model contains 27122 
equations, 818215 continuous and 1115 discrete variables. The CPU 
time spent to find a single Pareto solution ranges from 32.5 CPU s to 
30857.2 CPU s. The first one correspond to the extreme scenario 
that minimises the environmental impact for SoC equal to 27. The 
highest value is for the environmental impact minimisation 
scenario for SoC equal to 18. Generally, as the social factor increases, 
the calculation time decreases. In general, it is appreciated that the 
resolution of the economic criteria is less time consuming than the 
environmental criteria case. 


5.5. Results 

5.5.1. Objective Junctions optimisation 

Extreme scenarios are obtained by individual objective func¬ 
tions optimisation. They are characterised by the three set of values 
(NPV, Impact 2002+ and SoC) derived from NPV maximisation, 
Impact 2002 + minimisation and SoC maximisation. The optimum 
networks are represented by four types of matter flows that 
connect the different sites, disregarding the period (month of the 
year) when the flow appears in the resulting network. Those flows 
are: raw material from harvesting, stored raw material, dried 
matter and chipped matter, characterised by different colors in 
Fig. 5. Flowrates are obtained in tons/month. Note that any of the 
networks transports dried material. The chipped biomass is the 
most used alternative and the one that characterise all the 
networks. Since it has been assumed that cassava waste is produced 
into each community, proportional to the population, there is no 
need to employ an MV microgrid and to use the intermediate sites. 
Moreover, the high investment costs associated with an interme¬ 
diate site prevent their use. Notice that the simplest network comes 
out from the environmental impact minimisation, with only chip¬ 
ped matter. Table 4 lists the calculated capacity of the equipments 
installed at each site, for each objective function. 

Fig. 5(a) depicts the network that maximises the NPV. Some 
stored and non-stored raw matter are transported between loca¬ 
tions. The model proposes to install a generation plant in each 
community but allocates the needed pre-processing facilities in 
only four sites: Old Konkrompe, Fakwasi, Kumfia and Abamba (see 
Table 4). The four sites with pre-processing units are strategically 
located to best handle the biomass of all the communities. The 
minimum chipping capacity is installed (lOt/h) in almost all the 
communities, which is enough to process all the needed cassava 
waste. Therefore, Old Konkrompe, Fakwasi, Kumfia and Abamba are 
centralised processing sites of cassava waste, which are four of the 


Table 2 

Principal output values. 


Parameter Values 

r gasif (°C) 702.00 

Producer gas composition (on a mole basis) 

CO 23.93 

C0 2 10.49 

N 2 37.07 

H 2 20.88 

CH 4 3.58 

H 2 0 4.03 

Flowrate (kg/h) 35.33 

LHV (MJ/kg) 6.32 

CGE (%) 68.00 

Power (kW e ) 15.80 

i ) (%) 17.00 


five most populated communities. Fluxes of raw material exist to 
supply the largest communities. Pre-treatment units installed to 
minimise Impact 2002+ (Table 4) have its minimum capacity at all 
the sites, except at the two largest, Fakwasi and Kumfia, which 
relies on adapted dryers and G-ICE units. Old Konkrompe, Boniafo 
and Abamba also have gasification plants with capacities larger 
than the minimum. For SoC maximisation, again, some stored and 
non-stored raw matter are transported between locations (see 
Fig. 5(c)). Fluxes of raw material exist only to supply Bompa 
community. The model installs each type of unit at each location, 
but the capacities are different if compared with the previous case, 
obtaining as a result a more complex network. All the pre¬ 
treatment units’ capacities installed correspond to the minimum 
value, except for Kumfia, that has a larger dryer and collects chip¬ 
ped biomass from other communities. Again, the gasification units 
are scaled to meet their respective community demand. The results 
show that those networks require biomass storage in order to 
satisfy the constant electricity demand. The storage period ranges 
between 5 and 8 months. 

Table 5 summarises the three criteria evaluated for each 
optimal network. Note that the value of SoC is the same for the 
environmental friendly and for the most social favorable scenarios. 
The one that maximises the NPV when SoC is fixed to its 
maximum value has been identified as the maximum social 
criterion configuration. The economic value corresponding to the 
BSC configuration that maximises the NPV is $89895.95, which is 
the maximum NPV that can be obtained in this case study. This 
value is decreased by 59% in the environmental friendly network, 
and by 50% in the social network. For the environmental param¬ 
eter, it can be deduced that the difference among configurations 
mainly falls on transportation of chipped biomass. On the other 
hand, the Impact 2002 + value is reduced by 2.3%, if the LCA is 
optimised instead of the NPV. Note that if the most environmental 
friendly option is selected, the decision maker should be willing to 
compromise around 60% of the optimum NPV which implies just 
a 2% improvement in the environmental metric. In order to 
decrease the environmental impact, there is a necessity of more 
decentralised units. This requires a significant greater investment 
in comparison with the NPV optimum alternative. Notwith¬ 
standing, the gain in the environmental impact only results from 
a reduction in transportation. The NPV optimal solution has an 
environmental impact of 113.46 pts, whereas the Impact 
2002 + optimum scenario has 110.94 pts and the social optimum 
network has 111.43 pts. The environmental impact for the social 
BSC network is closer to the optimal environmental impact. The 
maximum value for the social criterion is 27, installing 3 units per 
site. The comparison between economic and social configurations 
show that: (i) in the former the capacity of the installed units is 
adapted to match the demand and, (ii) in the latter the units are 
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Table 3 

Economic parameters for pre-treatment units, G-ICE plant, transportation and utilities consumption. Data from [58,60,61 ], Ankur Scientific Energy Technologies PvtLtd., KITE and 
WACSO. 


Base scale Base investment 

MS10.5 
MSI.2 
M$0.05 a 

Tractor full $0.32/kmt 
$5000/km 

a LV network costs are included here. 
b Transformer cost is $1000. 


Drying 100 t/h 

Chipping 80 t/h 

G-ICE system 20 kW c 

Transportation biomass 
MV network b 


O&M (% of investment) Utility consumption Lifetime (y) 

3 0.06 • t Hl0ev (t diesel) 15 

20 Bond law 0.15 ■ t input (kW) 15 

4 7 
Loading and offloading $1.32/t 
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Fig. 5. Optimum networks configuration for the three selected criteria. Axes in km. See in brown the non-stored raw matter flow, in green, the raw matter flow, in orange the dry 
matter flow and in purple the chipped matter flow. 


installed and then operated to meet the demand. Thus, a certain 
centralisation is needed in DES to ensure the sustainability if the 
network. 

Fig. 6 depicts the investment and the annual costs breakdown 
for each described optimal BSC. Investment costs are $284422, 
$343187 and $331967 for NPV, Impact 2002+ and social optimised 


scenarios, respectively. The utilities costs for each BSC are around 
$26000/y due to the diesel consumption of the dryer, that processes 
about the same amount of biomass into each network. Note that the 
economic optimum scenario has the lowest investment costs, but 
the highest operational costs. The opposite occurs in the environ¬ 
mental friendly scenario. The difference results from the trade-off 


Equipment capacity f 


[ for the three selected criteria. 


Dryer (t/h) Chipper (t/h) G-ICE (MJ/h) Dryer (t/h) Chipper (t/h) G-ICE (MJ/h) Dryer (t/h) Chipper (t/h) G-ICE (MJ/h) 


Old Konkrompe 
Fakwasi 

Bompa 

Nwunwom 

Boniafo 

Abamba 























z-Fortes et at / Energy 44 (2012) 79-95 


Table 5 

Economic, environmental and social aspects for the individual objective 
functions optimised networks. 


NPV optimization ($) 89895.95 

Impact 2002+ (pts/y) 113.46 

Social criterion 17.00 


NPV ($) 36867.21 

Impact 2002 + optimisation (pts/y) 110.94 

Social criterion 27.00 


NPV ($) 45155.60 

Impact 2002+(pts/y) 111.43 

Social optimisation 27.00 


between transportation and decentralisation. The social optimum 
scenario is an intermediate network. 

Analogous observations can be retrieved from Fig. 7(a), where 
the environmental impact is depicted per SC echelon. The contri¬ 
bution that changes the most among the optimised networks is 
that of pre-processed biomass transportation. Biomass pre¬ 
processing and electricity generation constitute the major share 
among the intervention contributions; therefore the increase of the 
efficiency of these processes is crucial. Fig. 7(b) represents the 
environmental impact share per damage category, i.e. human 


health, ecosystem quality, climate change and resources. Those are 
in the same range for each optimum scenario. The highest impact 
comes from climate change, being a 80% of the whole impact. 

5.5.2. Pareto front 

The trade-off among the three selected parameters is studied in 
this section. SoC ranges between 17 and 27, and using the 
e-constraint 132 networks are found. After filtering, 104 networks 
result as non-dominated options. Fig. 8 representats the Pareto set 
of solutions, NPV vs. Impact 2002+, for each social iso-line. 

It is observed that as the SoC increases, the environmental 
impact decreases at the expense of compromising the NPV. The 
curves for SoC equal to 17 and 18 are far from the behaviour of the 
other level curves. Those are the better options with regard to only 
the economic criterion. In the considered SoC range, the NPV varies 
in an interval of M$53, while the environmental impact oscillates 
1.8 pts. They represent a decrease of 59% and an increase of 2% in 
comparison with the single optimisation values. 

5.5.3. Sensitivity analyses 

Several sensitivity analyses have been performed to show how 
optimum scenarios can vary if some key parameters are changed. 
The effect of the electricity price ($0.233/kWh) on the optimal NPV 


3 ■ Drying ■ Chipping ■ BG-GE ^ ■ Transport ■Utilities ■ O&M 



Investment share Annua! costs share 

Fig. 6. Breakdown of costs for NPV, impact 2002+ and social optimised networks. 



Fig. 7. Environmental impacts arising from single optimisation 


s, in impact 2002+ pts/y. 





















z-Fortes et at. / Energy 44 (2012) 79-95 



Impact 2002+ (pts) 

Fig. 8. Pareto set of solutions. 

network considering the maximisation of the NPV, is analysed. If 
the electricity price is decreased, then the viability of the network 
breaks down, i.e. the NPV starts to be below zero. Prices lower than 
around $0.2/kWh would require a subsidy policy. Fig. 9 shows the 
behaviour of the internal rate of return (IRR) when changing the 
electricity price for economic and environmental optimum 
scenarios. One can see that for the same price, the economic 
optimum scenario improves about 5% the IRR. For a IRR equal to 
zero, the maximum NPV configuration allows a electricity price 
value of $0.15/kWh; whereas, the most environmental configura¬ 
tion allows a price value of $0.16/kWh. 

5.5.4. Global remarks 

The needed amount of biomass to satisfy the electricity demand 
is 1526t/y. An 8.4% of the total available cassava rhizome is not used 
and can be employed for other purposes. The most important 
differences among the three selected criteria concern trans¬ 
portation and investment. Due to the cassava waste disposition, 
which is present in all the communities, no MV microgrid is 
installed. Therefore, the G-ICE plants capacity installed in all the 
communities, for all the Pareto solutions, are the same. The smallest 
communities whose demands are far from the biggest ones, such as 
Seneso or Nwunwom, are the communities that show more vari¬ 
ability along the different scenarios. The largest variability comes 
from the chipper and dryer installations: those are the units used to 
adjust the social factor during the search of optimum scenarios. 



networks. 


Notice that even though there is enough biomass to be self- 
sufficient, biomass from other communities is processed to supply 
the demand of small communities. By doing so, the investments in 
pre-treatment technologies are reduced. Such a reduction has 
a major contribution than the savings in transportation that would 
be gained if the biomass were pre-treated locally for these small 
communities. It is also relevant to notice that transportation of raw 
biomass is carried out from some small communities to other ones 
with installed pre-treatment capacity in order to face seasonality. At 
this point, the decision-maker is able to select one of the configu¬ 
rations represented in the PF, according to his/her criteria. 

6. Conclusions 

This paper has demonstrated the capability and adequacy of 
a multi-objective mathematical programming approach for the 
decision-making associated with a DES network design. The 
mathematical model comprises three performance metrics with 
economic, environmental and social concerns. There does not exist 
a unique modelling and global approach to tackle such problems: 
the strategies must be adapted depending on the raw material, 
disposable pre-treatment units and final energy purposes. 
Distributed approaches should rely on a certain level of central¬ 
isation to be feasible along the planning horizon and favour the 
economic and social developments of regional areas. The capabil¬ 
ities of the proposed model have been exemplified with a case 
study placed in Ghana. This case study showed some centralisation 
of pre-treatment equipment to take economic advantages that may 
help to the project feasibility. The multi-objective optimisation 
resulted in a PF that provides crucial information for the design 
problem. It showed the main trends of the optimal solutions, 
revealing trade-offs among the economic criterion, environmental 
and social criteria. Once the model is built, it is possible to perform 
sensitivity analyses to provide information about sensible param¬ 
eters and their influence in the network design. In this paper, the 
electricity price is analised, being eventually useful to establish 
price and subsidies policies that can help to motivate such projects. 
All the compiled information is an important input for the decision¬ 
making process and support the selection of the most appropriate 
supply chain configuration. 

The energy block is the most versatile part of the model in DES 
configuration: this part can be easily adapted to other micro¬ 
generation sources, such as wind, hydro or solar, or even to 
investigate the use of more than one renewable source in hybrid 
systems. Part of our research agenda includes the evaluation of co¬ 
generation and tri-generation to enhance the supply chain effi¬ 
ciency. The present approach can be enhanced to consider other 
advanced pre-treatment methods, such as pelletisation, torre- 
faction or pyrolysis. The mathematical model presented in this 
work deals with the combination of feeds and determines the 
optimum proportions of each input for the energy generation 
activity. Further work is devoted to the extension of this flexibility 
to other activities included in the biomass supply chain, such as 
storage and pre-treatment. Decomposition strategies to reduce the 
computational cost of solving this type of problems is also devised 
as future work. 
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Notation 


Superscripts and subindices 
ar as received basis 

e electric 

th thermal 

ev evaporated 

Acronyms 

BD bulk density 

BM biomass 

BSC bio-based supply chain 

CGE cold gas efficiency 

COE cost of the energy 

d day 

DES distributed energy systems 

DM dry matter 

ER equivalence ratio 

G gasification 

ICE internal combustion engine 

1RR internal rate of return 

LCA life cycle assessment 

LCI life cycle inventory 

LCIA life cycle impact assessment 

LHV lower heating value 

LP linear programming 

LV low voltage 

MC moisture content 

MFP multi-functional platform 

MO-MILPMulti-objective mixed integer linear program 

MV medium voltage 

NPV net present value 

O&M operation and maintenance 

PF Pareto front 

SC supply chain 

SCM supply chain management 

STN state task network 

Indices 

a mid-point environmental impact categories 

e suppliers 

/,/' facility locations 

g end point environmental impact categories 

i tasks 

j equipment technology 

k intervals for piecewise linear approximation of 

economies of scale 
s materials (states) 

t, t' planning periods 

Sets 

A g set of mid-point environmental interventions that are 
combined into endpoint damage factors g 
Em, set of suppliers e that provide raw materials 

Eprod set °f suppliers e that provide production services 

Etr set of suppliers e that provide transportation services 

FP set of materials s that are final products 

/ set of tasks i with variable input 

I] set of tasks i that can be performed in technology j 


J e technology j that is available at supplier e 
]f technology j that can be installed at location/ 

Ji technologies that can perform task i 

Jstor technologies to perform storage activities 
Mkt set of market locations 

NTr set of production, or non-transport, tasks 
RM set of materials s that are raw materials 
Sstor set of materials/biomass that if stored change their 
properties 

Sup set of supplier locations 

T s set of tasks producing material s 

T s set of tasks consuming material s 

Tr set of distribution tasks 

Parameters 

Asft [kg] maximum availability of raw material s in period t in 
location/ 

Dem s f t [MJ] demand of product s at market/in period t 
distance•%, [km] distance from location /to location/' 

FCFJjji [$/h] fixed cost per unit of technology j capacity at location/ 
in period t 

FEjjfr' 1 [ h ] increment of capacity equal to the upper limit of interval k 
for technology j in facility/ 
i r [adim.] discount rate 

Invest [$/km] investment required for medium voltage 
connectivity per length unit 
M [adim.] a big number 

NormFg [adim.] normalising factor of damage category g 
Pricey [$/MJ] price of product s at market /in period t 
Price^ l,t [$/MJ] investment required for an increment of capacity 
equal to the upper limit of interval k for technology 
j in facility / 

Tortuosity [adim.] tortuosity factor 
Water s [adim.] Moisture for material s 
Water? ,ax [adim.] Maximum moisture for task i performed in 
equipment j 

Greek symbols 

a S ij [adim.] mass fraction of task i for production of material s in 
equipment j 

a s y [adim.] mass fraction of task i for consumption of material s in 
equipment j 

py [adim.] minimum utilisation rate of technology j capacity that is 
allowed at location / 

!', ag [adim.] g endpoint damage characterisation factor for 
environmental intervention a 

%. [h/kg] capacity utilization rate of technology j by task i whose 
origin is location/and destination location/' 

Ptfgt [$/kg] unitary transportation costs from location/to location 
/' during period t 

z ijfet [S/kg] unitary cost associated with task i performed in 

equipment j from location/and payable to external 
supplier e during period t 

T sfet [S/kg] unitary cost associated with handling the inventory of 
material s in location/and payable to external supplier e 
during period t 

Xest [$/kg] unitary cost of raw material s offered by external 
supplier e in period t 

i'ijffa [points/kg] a environmental category impact CF for task i 

performed using technology j receiving materials 
from node/and delivering it at node/' 
yp]j a [points/(kgkm)] a environmental category impact CF for the 
transportation of a mass unit of material over 
a length unit 



Binary variables 

Vjft 1 if technology j is installed at location/in period t, 

0 otherwise 

Zff, 1 if facilities f and f are interconnected by a medium 
voltage line, 0 otherwise 

SOS2 variables 

Zjfkt variable to model the economies of scale for technology j 
in facility/at period t as a piecewise linear function 


Continuous variables 

DamCgft [point] normalised endpoint damage g for location/in 
period t 

DamCg C [point] normalised endpoint damage g along the whole SC 
EPurchet [$] economic value of purchases executed in period t to 
supplier e 

ESalest [$] economic value of sales executed in period t 
FAsset t [$] investment on fixed assets in period t 
FCostt [$] fixed cost in period t 

Fjft [h] total capacity of technology j during period t at location/ 
period t 

IC a ft [point] mid-point a environmental impact associated to site/ 
which rises from activities in period t 
Impactj 002 [point] total environmental impact for site/ 
Impact™^ [point] t 
NPV [$] econom 

Pijff't [kg] specific activity of task i, by using technology j durir 
location/' 

Profitt [$] profit achieved in period t 

Pvsijft [kg] input/output material of material s for activity of t£ 

period tin location/(This n 
Profitt [$] profit achieved in period t 
Purch p / t [$] amount of money payable 
associ, 

Purch^p [$] amor 


Purch" [$] amount of 
associated 
Sales sS , t [$] amount of 
period t 

Ssft [kg] amount of sti 
SoC [unit] surrogate sc 



e & Commercial. It refers tc 
of electricity, i.e. motors, large fridges or freezers, agro¬ 
processing or sawmill equipment. 

To estimate the electricity demand in Atebubu district, we have 
considered a reference community. It is a small rural community 
without previous access to electricity with 10 potential users and 
with 7 households. For further details on values please see [65], 
Using the reference demand, the net demand of each commu¬ 
nity from our case study has been estii 
number of households. The current a 
will be supplied with a distribution n 
hence, the intra-community L.V distril 
mated to approach the gross demand, that is, the electricity that 
needs to be supplied to each community. Based on a related 
International Electrotechnical Comission (1EC) standard [66], the 
maximum losses in an LV monophasic distribution microgrid 
would be kept at 5%, plus a 1% of self-consumption due to the 
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Appendix A. Data for tl 

Three main types of 
follow the typical regulat 
reference ii 
the rural electrii 



• Residential. It 
as lighting, 
communities 
households (up to 70% of 
tricity demand, in the range of 8 to 20 kWh [ 
(equivalent to 250 to 700 Wh/d), even with access 
service scheme. 

• Community. It refers to the use of electricity in con 
public premisi 
water pumpin 
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